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We evaluate the upward shower and muon event rates for two characteristic four neutrino mixing 
models for extragalactic neutrinos, as well as for the atmospheric neutrinos, with energy thresholds 
of 1 TeV, 10 TeV and 100 TeV. We show that by comparing the shower to muon event rates, one 
can distinguish between oscillation and no-oscillation models. By measuring shower and muon event 
rates for energy thresholds of 10 TeV and 100 TeV, and by considering their ratio, it is possible 
to use extragalactic neutrino sources to determine the type of four-flavor mixing pattern. We find 
that over several years of data taking, a kilometer-size detector has a very good chance of providing 
valuable information about the physics beyond the Standard Model. 



I. INTRODUCTION 

The combined results of solar, atmospheric and labo- 
ratory experiments with neutrinos, taken at face value, 
require a fourth neutrino species. This follows from the 
observation that the results of the three categories of 
experiments require at least three mass-squared differ- 
ences 5m^. The mass-squared difference for solar neu- 
trino experiments is limited to Stti^^i^,^ ^ 10""^ eV^ [Q. 
The SuperKamiokande results ||] for atmospheric neutri- 
nos require Sm^f.^ ~ 3 x 10^^ eV^, and the laboratory 
LSND experiment |] limits 5m\s^jj > 0.2 eV^. While 
precision measurements of the invisible decay width of 
the Z'^ boson constrain the addition of a fourth genera- 
tion neutrino species with weak interactions, they do not 
constrain the possibility of a sterile neutrino species mix- 
ing with the ordinary electron, muon and tau neutrinos. 
Analyses of solar and atmospheric experiments assum- 
ing pure Vi Vs oscillations indicate that the data do 
not support the hypothesis, however, combined fits to 
all of the experimental data do require a sterile neutrino 
species [Q. 

Recent SuperK data on atmospheric neutrinos indicate 
Vu ^ Vt oscillations with mixingbeing nearly bi-maximal 
[01 . We have shown in Ref. [g| that for extragalactic 
sources of muon neutrinos, mixing with tau neutrinos in 
transit to the Earth leads to distinct signatures of oscil- 
lation which do not require explicit identification of a tau 
lepton in large underground experiments. In this paper, 
we investigate the signatures of neutrino oscillations in 
large underground experiments for models with three ac- 
tive neutrinos and one sterile neutrino whose mixing is 
constrained by lower energy experiments. 



with (5m^t„ and a splitting between the two nearly degen- 
erate pairs characterized by 5'm\gj^jj. The approximate 
mixing matrix in this scenario has the form 
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where e < 0.1 Following Ref. we caU this a 2 -I- 2 
scenario. 

The second pattern of neutrino mass and mixing, al- 
lowed by the most recent analysis of the LSND experi- 
ment, has large mixing between the three ordinary neu- 
trinos and a small mixing with the sterile neutrino, char- 
acterized by small parameters e and 8. A characteristic 
mixing matrix has the form 
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Following Ref. jj], we call this a 1 -I- 3 scenario. 

For extragalactic sources, because of the large dis- 
tances involved, the oscillation flavor ratios are essen- 
tially independent of mass-squared differences and neu- 
trino energy, since 
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II. MIXING MODELS AND EXTRAGALACTIC 
FLUXES 

Global fits to oscillation data fall into two distinct pat- 
terns of neutrino mixing. The first is a mass spectrum 
in which i/^ ^ Vg with a mass splitting Sml^i^^, ^ 



for 8m} in eV^, L in km, and E in GeV. Consequently, 
the oscillation probabilities can be written as 
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in terms of the elements of the neutrino mixing matrix 
Uij. The i/g flux in detectors F[' , in terms of the fluxes 
at the source Ff are [^,0] 
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As we discuss below, our "source" fluxes are actually 
summed over many sources to yield isotropic fluxes. We 
continue to denote these isotropic fluxes, unmodified by 
oscillations, as the source fluxes. 

Given source ratios of fluxes : : : ^ Q : 
1:2:0 yield different ratios of neutrino fluxes at the 
detector, depending on whether the 2 + 2 or 1 + 3 scenario 
describes four-neutrino mixing. In the 1 + 3 case with 
small e and 5, the detector ratios are approximately : 
1:1:1, while for the 2 + 2 case, sterile neutrinos make 
an important component of the flux at the Earth, with 
r/f : r/f : r/^ : i/f ~ 0.5 : 0.5 : 1 : 1. With three neutrino 
species and bi-maximal mixing, one flnds : : ~ 
1:1:1. Without mixing, the source fluxes and detector 
fluxes have the same flavor ratios of 1 : 2 : 0. As a result, 
the 1 + 3 scenarios essentially reproduce the three-flavor 
bi-maximal mixing scenario, while the 2 + 2 scenarios lie 
between the 3-flavor bi-maximal mixing model and the 
no- mixing model. 

There are a variety of predictions for isotropic neu- 
trino fluxes from active galactic nuclei (AGN), gamma 
ray bursters (GRB) and models with topological defects. 
A sample of these predictions for muon neutrino plus an- 
tineutrino fluxes, in the absence of oscillations, arc shown 
in Fig. 1. These include the AGN models of Stecker 
and Salamon (AGN_SS) g and of Mannheim Model A 
(AGN_M95) 0. Both of these models predict neutrinos 
fluxes that represent the upper bounds for their class of 
the models. In particular, the Stecker-Salamon flux is an 
upper bound for AGN core emission, while Mannheim 
Model A is an upper bound for AGN jet emission mod- 
els. The Stecker-Salamon flux is bound by the the dif- 
fuse X-ray background, while Mannheim flux is bound 
by the extragalactic gamma ray background. GRB pre- 
dictions are represented by the Waxman-Bahcall model 
of Ref. (GRB.WB). Two topological defect models, 
from Sigl, Lee, Schramm and Coppi (TD.SLSC) (l|] and 
Wichoski, MacGibbon and Brandenberger (TD_WMB) 
Ip^ are also shown. We include fluxes with energy be- 
haviors like 1/E and 1 /E"^ as well. We chose normaliza- 
tion for each of these fluxes consistent with the current 
limits [0, F^^ ^ 10-12 (£;/GeV)-iGeV-icm-2sr-is-i 

and = 10-*^ (i;/GeV)-2GeV-icm-2sr-is-i (in 

which tfie fluxes include neutrinos and antineutrinos in 
equal amounts and correspond to the initial fluxes be- 
fore accounting for oscillations). Our choice of normal- 
ization for E^^ flux is a factor of ten larger than re- 
cently proposed neutrino flux associated with GRB ( 
F^' = 10-^ (£;/GeV)-2GeV-icm-2sr-is-i) (w), and 
a factor of 50 larger than the upper bound for strong 
source evolution previously discussed by Waxman and 
Bahcall |l^]. The E~^ flux is smoothly cut off at high 
energies, as indicated in Fig. 1. 

Although the normalizations of the fluxes of neutri- 
nos from astrophysical sources are uncertain, we evaluate 



event rates quantitatively for these representative energy 
behaviors of the fluxes which, as models improve, can be 
rescaled to accommodate different normalizations. 
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FIG. 1. Isotropic muon neutrino plus antineutrino flux pre- 
dictions for AGN models (solid lines, upper curve at low en- 
ergy corresponds to AGNJVI95, while the lower curve is for 
AGN_SS model) , GRB (dotted line) , topological defects mod- 
els (dash-dotted Unes, upper curve corresponds to TD_(Model 
A) WMB, while the lower curve is for TD_SLSC), E'^ flux 
(lower dashed line at low energy) and (upper dashed line 
at low energy) and angle-dependent atmospheric (ATM). The 
flux is scaled by neutrino energy squared and the antineutrino 
flux is taken equal to the neutrino flux. 

For comparison, the hatched curve shows the angle de- 
pendent atmospheric flux (ATM) ||l^ from kaon and pion 
decays (conventional flux). We have extrapolated this 
flux beyond the 10 TeV range given in Ref. ||l^ using 
an angular dependent power law. This does not account 
for the change in the input cosmic ray spectrum, the so- 
called "knee" at energies of - 10^ GeV 11, which trans- 
lates to a reduced conventional neutrino flux at high en- 
ergies. Thus, we overestimate the conventional neutrino 
event rates, especially at our highest threshold of 100 
TeV, however, at that threshold, the atmospheric rates 
are small. Our evaluation does not include the prompt 
neutrino flux from semileptonic decays of charmed parti- 
cles produced in the atmosphere, which result in neutrino 
fluxes with a power law increased by one factor of energy. 
Recent evaluations of the prompt neutrino flux suggest 
that it is important only above 100 TeV |p^|-pOt, although 
it has been emphasized that there are large theoretical 
uncertainties in the evaluation |2^ . 

In our previous work we have evaluated how neu- 
trino interactions in the Earth modify these representa- 
tive neutrino fluxes with the regeneration attributed to 
the neutral current interactions. In addition, for tau neu- 
trinos, we described the extent to which neutral current 
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interactions and charged current production of r followed 
by its decay, regenerate neutrinos. There, we made a de- 
tailed numerical evaluation of the "pile-up" of tau neu- 
trinos that was discussed by Halzen and Saltzberg in Ref. 
p2t . In the next section, we use the neutrino fluxes of 
Fig. 1 modified due to their passage through the Earth 
in the appropriate flavor proportions for the 1 + 3 and 
2-1-2 scenarios. 



III. UNDERGROUND SIGNATURES 

Backgrounds from atmospheric muons make downward 
event rates difficult to extract, so our focus is on upward 
events. Two types of events will be produced: muon 
events and shower events. The muon events come from 
upward muons from ^ ^ charged current events and 
from charged current production of taus followed by a 
muonic decay, v^. ^ t ^ fiX. In spite of the pile-up in 
the tau neutrino flux, the net effect of oscillations is to 
reduce the muonic event rate by approximately a factor 
of two relative to the no-oscillation rate, whether in the 
1 + 3 or 2 + 2 mixing scenario. Given the uncertainties in 
the normalizations of the extragalactic fluxes (unlike the 
atmospheric flux), the muon event rate is not enough to 
distinguish oscillation scenarios from no-oscillation sce- 
narios. 

The distinction between oscillation and no-oscillation 
scenarios comes by comparing the muon rate with the 
shower rate. The principle of this procedure is very simi- 
lar to what some experiments will be able to do with low 
energy neutrinos: namely, to compare the neutral cur- 
rent rate to the charge current interaction rate. In large 
neutrino telescopes, tau identification is difficult, and we 
assume that showers of hadronic and electromagnetic ori- 
gin cannot be distinguished. Consequently, the shower 
rate comes from the following processes from neutrino- 
nucleon (N) interactions: 

Vt-N —t T + hadrons, r — > i/,- -|- hadrons , 
Vt-N — > t + hadrons, t ^ Vt + e + Vg , 

Vt-N I'r + hadrons, (6) 
Ve,p,N Ve^fi + hadrons, 
VgN ^ e + hadrons . 

The shower energy is taken as the sum of hadronic and 
electron energies including both the production vertex 
and the decay vertex, where applicable. 

The upward muon rates are determined from an eval- 
uation of the quantity 

Rate^ ANaJ^ dE, J dy j dz{R^(E^,Ef''))'^ (7) 

X (E,, ,X)e{E^- Ef-) {E,,y, z) 

where i?^ = Ei,{l — y)z. The function / depends whether 
the source of the muons is muon neutrino charged current 



interactions or tau neutrino charged current interactions 
followed by the muonic decay of the tau: 

ff,{E„,y, z)= d{l - z) ioi- i^f, ^ , (8) 

ff,{E^,y, z)= — for i/^ ^ T ^ pi , 

where the decay formula appears in the appendix of Ref. 

. The quantity A is the effective area (taken as 1 km^ 
here), Na is Avogadro's number and the differential cross 
section is for neutrino-interactions. The average range of 
the muon {R) for initial energy E^ and final energy 
i?™'" is the additional space dimension that makes up 
the target volume. The column depth {X) dependence 
of the neutrino flux represents the angular dependence of 
the flux due to the attenuation of the flux after passage 
through the Earth, and in the case of atmospheric neu- 
trinos, includes the angular dependence of the flux at the 
surface of the Earth. 

For the shower rate, one does not have the benefit of 
the muon range, so the full instrumented volume V enters 
into the calculation of the event rate, here taken to be 1 
km^: 

Rate= VNa J dE^ J dy J dz^ (9) 

shr 

xF,(i?„ X)e(£;,hr - SX") UE,, y, z) , 

where fs has a similar form to . Details of the differen- 
tial decay distribution of the tau used in the evaluation, 
as well as the specific dependence of i?shr on E^, y and 
z appear in The differential cross section is either 
charged current or neutral current, depending on which 
process in Eq. (6) is being considered. 

The neutrino and antineutrino cross sections are 
evaluated using the CTEQ5 parton distribution functions 
pBI , and the attenuation of the fluxes assume Earth den- 
sities of the Preliminary Earth model described in Ref. 
p6| . As in Ref. our upper bound of integration is 10^ 
GeV because the attenuation of Vt- fluxes at higher en- 
ergies have not been evaluated, due to the complication 
of including tau energy loss at high energies ||2^. Our 
results for the 100 TeV threshold are therefore conserva- 
tive. 

The ratio of the shower rate to the muon rate for the 
extragalactic flux models of Fig. 1 are used to bracket 
the theoretical expectations for the ratio in 2-1-2 and 1-1-3 
mixing models, and compared with the standard model 
expectations. As in Ref. we separate the flux mod- 
els into two categories: more steeply falling fluxes (i?~^, 
ATM and the Mannheim AGN flux) and the less steep 
behavior of the other five sample fiuxes of Fig. 1. Exper- 
imentally, the steeply falling fluxes and the models with 
a weaker energy dependence are characterized by distinct 
muon event rate nadir angle dependence and the thresh- 
old energy dependence. The two separate categories of 
fiuxes, for threshold energies of 1, 10 and 100 TeV are 
shown in Figs. 2 (a-f). The 1 + 3 models are bracketed 
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by the black lines, the 2 + 2 models are represented by 
the shaded area and the standard model with the dark 
band. 

For the 1 TeV threshold, the 1 + 3 and 2 + 2 event rates 
as a function of nadir angle overlap for both energy be- 
haviors of the input fluxes. The no-oscillation band lies 
below the oscillation bands with the exception of the near 
horizontal rates in the steeply falling flux category (la- 
beled E^"^). The separation between the oscillation and 
no-oscillation scenarios is better for higher thresholds, 
looking only in terms of the theoretical ratio of shower 
to muon event rates. However, the higher thresholds have 
lower event rates overall. 

The bands in Fig. 2 indicate the different ratios for 
oscillation models relative to no-oscillation models, how- 
ever, statistical errors from the signal as well as from the 
atmospheric background need to be included in order to 
interpret whether or not oscillation scenarios are distin- 
guishable. In the next section, we include a discussion of 
statistical errors. 



IV. DISCUSSION 

In Tables Till we present our results for the shower 
and muon event rates integrated over the nadir angle for 
energy thresholds of 1 TeV, 10 TeV and 100 TeV and 
for a kilometer-size detector. The no-oscillation results 
agree with evaluations in Ref. ]2^ for fluxes in common, 
modulo corrections due to different parton distribution 
functions, a more precise evaluation of neutrino attenua- 
tion in the Earth ||2^ and the upper limit of energy inte- 
gration. We do not include the topological defect models 
in these tables because the event rates are so low. For the 
1 TeV threshold, only 0.1 — 0.2 upward muon events are 
predicted per year for a km^ instrumented area for the 
TD_WMB model. In all cases, we have considered shower 
rates in a km'^ volume, and muon rates/km^. The two 
oscillation scenarios give the same muon event rates but 
very different shower rates. From Table I, we note that 
for 1 TeV energy threshold, atmospheric background is 
large. 
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FIG. 2. Ratio of shower event rates to muon event rates 
for the 1 + 3, 2 + 2 and no oscillation scenarios as a function 
of nadir angles for threshold energies of 1 TeV, 10 TeV and 
100 TeV, for the indicated fluxes. The black lines bracket the 
1 + 3 scenario, the shaded area represents the 2 + 2 scenario 
and the dark band represents no-mixing case for the given 
fluxes. 
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FIG. 3. Integrated over solid angle, the ratio of upward 
shower to muon event rates for given flux models for a thresh- 
old energy of 1 TeV. The error bars for 1 + 3 scenario (top 
data), 2 + 2 scenario (middle data), and no-mixing case (bot- 
tom data) are determined assuming n years of data taking. 

In Figs. 3-5, we show the ratio R of the shower to muon 
event rates integrated over nadir angles {R = Nshr/N^) 
for energy thresholds of 1 TeV, 10 TeV and 100 TeV for 
a variety of fluxes and for different oscillation scenarios. 
The ratio of event rates in Figs. 3-5 have the atmospheric 
flux subtracted. The errors in the shower and muon rates 
(ashr and tr^) are evaluated with the assumption of a 
Poisson distribution, including the statistical error on the 
number of background events from the atmospheric flux. 
The error in the ratio of event rates is then given by. 
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(TR = RJ — h 



shr 



N„, 



(10) 



The event rates Nshr and are from Tables 1-111, scaled 
by n years as indicated in Figs. 3-5. 

From Fig. 3 and Tabic 1, wc note tliat AGN fluxes, as 
well as E~'^ and fluxes, would have the best possibil- 
ity of separating different oscillation scenarios by detect- 
ing the showers and muons in kilometer-size detector over 
the period of one to six years, assuming a good under- 
standing of the atmospheric muon neutrino and electron 
neutrino fluxes and normalizations near their current up- 
per bounds. After one to six years of data taking, one 
would be able to separate oscillation from no oscillation 
scenario for the fluxes in Fig. 3. The GRB_WB flux 
would require more than twenty years of data to reduce 
statistical errors so that the error bars do not overlap for 
the two oscillation scenarios. 



1.5 



B 1.0 

ca 
K 

B 

> 



0.5 - 



0.0 







I 




10 TeV 




E"' 


I 


AGN_SS 






I 




I 
i 






I 










i 




ATM 






n = 1 








i 










n = 3 




n = 1 


n = 1 



FIG. 4. Integrated over solid angle, the ratio of upward 
shower to muon event rates for given flux models for a thresh- 
old energy of 10 TeV.The error bars for 1 -f 3 scenario (top 
data), 2-1-2 scenario (middle data), and no- mixing case (bot- 
tom data) are determined assuming n years of data taking. 

The ratio of shower to muon event rates for an energy 
threshold of 10 TeV are shown in Fig. 4. As seen in Table 
II, the AGN model of Stecker and Salamon predicts large 
shower and muon rates, an order of magnitude larger 
than the atmospheric background. The AGN model of 
Mannheim and the GRB model of Waxman and Bah- 
call predict rates to be about 15-40% of the background, 
while E~'^ and fluxes give significantly more events 
than the atmospheric background when oscillations arc 
taken into account, assuming that the fluxes occur with 
a normalization at their current upper limit. Thus, there 
is a possibility of detecting extragalactic neutrinos by im- 
posing a 10 TeV energy threshold [5]. We note that for 
the AGN_SS flux, as weU as for E~'^, E'"^ fluxes, one 



to three years would be sufficient to separate different 
oscillation scenarios. For the AGN_M95 and GRB.WB 
fluxes, it would be necessary to take data for seven to 
nine years respectively. The 10 TeV threshold is gen- 
erally more favorable than the 1 TeV threshold because 
of the reduction of the atmospheric background and a 
sufficiently large signal. 

If the energy threshold is increased to 100 TeV, the at- 
mospheric background is small. The shower and muon 
event rates are several hundreds for AGN_SS model, 
while for all other models the rates are much smaller, 
but still significant, for kilometer-size detector in a time 
period of one year. Furthermore, the ratio of the shower 
and muon events can clearly separate oscillation and no- 
oscillation scenarios, as well as distinguish two oscilla- 
tions scenarios that we consider. In Fig. 5 we show this 
ratio including statistical errors. Wc find that even with 
error bars it still seems possible to distinguish between 
different oscillation scenarios. However, in order to re- 
duce statistical errors, the GRB.WB and AGN_M95 flux 
predictions would require collecting data for more than 
a decade. On the other hand, one year is sufficient for 
the AGN_SS flux. In the case of an E~^ flux, we expect 
several hundreds of showers and muons over the period 
of three years, which is necessary in order to reduce error 
bars to the point of being able to distinguish oscillation 
scenarios. For a steeper flux, E~'^, one year is sufficient 
to get enough events as well as to reduce the statistical 
errors at the current AMANDA limit on the isotropic 
flux normalization. 
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FIG. 5. Integrated over solid angle, the ratio of upward 
shower to muon event rates for given flux models for a thresh- 
old energy of 100 TeV. The error bars for 1-1-3 scenario (top 
data), 2-1-2 scenario (middle data), and no-mixing case (bot- 
tom data) are determined assuming n years of data taking. 

Figs. 3-5 show the feature illustrated by Fig. 2, when 
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specific fluxes are integrated over nadir angle. Fluxes 
with similar spectra cluster with the same ratio, here, 
the Stecker-Salamon flux and E^^ have similar energy 
behaviors and similar ratios. Therefore, by using muon 
event energy and nadir angle dependence to characterize 
the energy spectrum of the incident flux, even without a 
detailed knowledge of the source of neutrinos, the ratio 
R will point to whether or not muon neutrinos mix with 
tau neutrinos in one of the two schemes described here. 

The ratios in Figs. 3-5 differ for the different energy 
dependences of the fluxes for several reasons. The muon 
range increases with energy, so the muon event rates are 
sensitive to the energy dependence in a different way than 
the shower rates which have a constant volume factor. 
For less steep fluxes, e.g., E~-^ and AGN_SS, the high 
energy part of the spectrum gives a larger contribution 
to the muon rate than for fluxes with a steeper energy 
behavior, e.g., . The ratio R of shower to muon 
event rates tends to increase with the steeper fluxes, and 
tends to increase with energy threshold. For all except 
the atmospheric flux, the electron neutrino contribution 
to the shower rate is substantial. The electron neutrino 
component of the atmospheric flux is small because of the 
high energies considered here, so the atmospheric ratio is 
not directly comparable to the ratios for other fluxes. 



V. CONCLUSION 

We have shown that measurements of the upward 
showers and muons, and considering their ratio, with 
the kilometer-size detector, could provide very good 
test of different oscillation scenarios. Taking into ac- 
count statistical errors, we find that generic fluxes, such 
as F^^ = 10-12 (£;/GeV)-iGeV-icm-2sr-is~i and 

F^^ = 10-^ (i;/GeV)-2GeV-icm-2sr-is-i, with en- 
ergy threshold of 10 TeV and 100 TeV, would have 
enough statistics in one to four years of data taking to 
distinguish not only between oscillation and no oscillation 
scenario, but also whether it is 2-f 2 or 1-1-3 model. The 
AGN model of Stecker-Salamon with a 10 TeV energy 
threshold predicts between 1200 and 2000 muons/km^ 
per year and between 500 and 800 shower events/km'^ 
per year, with clear separation between different oscil- 
lation scenarios. The other AGN model, proposed by 
Mannheim predicts less events and would require 
seven years to reduce statistical errors in the ratio of 
shower to muon event rates such that there is clear dis- 
tinction between not only the no oscillation and oscil- 
lation cases, but also different oscillation models. If we 
want to use GRB_WB flux to test these models, it 
would take nine years with the kilometer-size detector to 
get sufficient statistics. A higher energy threshold of 100 
TeV makes the atmospheric neutrino background small, 
but event rates for extragalactic neutrinos are also re- 
duced. Still, the AGN_SS model would be able to pro- 
vide valuable information about the oscillation scenario, 



as well as E flux in just one year of data, while an -"^ 
flux with our normalization would need four years. Other 
ffuxes, such as GRB_WB and AGN_M95 flux would re- 
quire much longer time. 

The assumption of the existence of sterile neutrino is 
necessary to explain the solar, atmospheric, reactor and 
accelerator data. The MiniBooNE experiment will 
search for j/^ — > oscillations and test the LSND results. 
The presence of the sterile neutrino in the 2-1-2 model 
can be tested by the measurement of the suppression of 
NG/CC ratio of the solar neutrinos by the SNO experi- 
ment whereas 1-1-3 model can be tested by searches 
from small amplitude oscillations at short baseline exper- 
iment such as ORLaND Since the 1-1-3 scenario at 
high neutrino energy mimics the three flavor bi-maximal 
mixing between <-> v^, the results here can also be 
applied to three-flavor models. 

Observation of the flavor ratio of extragalactic neutri- 
nos could serve as a complementary test for oscillation 
models that incorporate sterile neutrinos. Oscillations 
of neutrinos from astronomical sources are averaged over 
very long baseline and thus cannot provide any direct in- 
formation about the neutrino mass. However, we have 
shown that extragalactic neutrinos could be used to dis- 
tinguish between different oscillation scenarios. The ob- 
served muon and shower event distributions as a function 
of angle and energy threshold will help determine the 
energy behavior of the incident extragalactic neutrino 
flux and point to a model or class of models for their 
sources. Combined measurements of the upward showers 
and muons and, in particular, their ratio provides a basic 
test of the oscillation scenario. 
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Note added: 

After the submission of our paper, the SNO collabo- 
ration has reported their first results [|32| which together 
with the SuperKamiokande results [|1| favor oscillation 
into active flavors, however, given the theoretical and ex- 
perimental uncertainties, oscillations into sterile neutri- 
nos are still a possibility. 
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TABLE I. Integrated upward shower (muon) rates per yeax for energy threshold of 1 TeV. 


Model 
1+3 

2+2 
No osc 


170(422) 
141(422) 
88.1(646) 


E-'^ AGN_SS AGN3195 GRB.WB ATM 
2528(3255) 918(2521) 224(157) 26.9(65.8) 
1898(3255) 726(2521) 163(157) 20.9(65.8) 
1750(5676) 593(4228) 158(280) 18.3(113) 2355(7346) 


TABLE II. Integrated upward shower (muon) rates per year for energy threshold of 10 TeV. 


Model 
1+3 
2+2 
No osc 


150(217) 
124(217) 
75.4(342) 


E-'^ AGN_SS AGN^95 GRB.WB ATM 
954(875) 809(1245) 34.4(20.2) 21.3(28.4) 
706(875) 633(1245) 24.5(20.2) 16.3(28.4) 
652(1540) 514(2139) 24.4(36.3) 14.2(49.7) 58.5(246) 




TABLE III. 


Integrated upward shower (muon) rates per year for energy threshold of 100 TeV. 


Model 
1+3 

2+2 
No osc 


84.8(54.2) 

68.4(54.2) 
42.0(91.3) 


E-'^ AGN.SS AGNJvI95 GRB.WB ATM 
182(100) 386(255) 2.46(0.999) 7.28(4.03) 
130(100) 285(255) 1.69(0.999) 5.21(4.03) 
124(182) 251(460) 1.78(1.83) 4.99(7.31) 0.774(4.18) 
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